Scaling and Universality in the Anisotropic Kondo Model and the Dissipative 

Two— State System 
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Scaling and universality in the Ohmic two-state system is investigated by exploiting the equiv- 
alence of this model to the anisotropic Kondo model. For the Ohmic two-state system, we find 
universal scaling functions for the specific heat, Ca(T), static susceptibility, Xa{T), and spin relax- 
ation function Sa(ijj) depending on the reduced temperature T / Ar (frequency tj/Ar), with Ar the 
renormalized tunneling frequency, and uniquely specified by the dissipation strength a (0 < a < 1). 
The scaling functions can be used to extract a and A,- in experimental realizations. 
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The low energy effective model for a large number of 
physical systems corresponds to a two-state system cou- 
pled to an environment . Examples include two-level 
atoms coupled to the electromagnetic field in quantum 
optics, electron-transfer reactions in biological systems, 
and the tunneling of defects in metallic glasses The 
simplest model for describing such systems is the spin- 
boson Hamiltonian [Q, 

HsB = -^Ao-j, + ieo-^ + ^ Wi(a|ai -I- 
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Here, the Pauli matrices ct^, i — x,y,z describe the two- 
level system, A is the bare tunneling matrix element be- 
tween the states tr^ —] and Uz —[ and e is a bias. The 
environment is represented by an infinite set of harmonic 
oscillators (labeled by the index i) with masses rui and 
frequency spectrum uji linearly coupled to the coordinate 
Q = \(lo^z with coupling constants Ci. In this paper we 
restrict ourselves to Ohmic dissipation, for which the en- 

vironment spectral function J(w) — § X]i( in L- 
takes the form J{uj) = 27rau;, for lo ^ a;,,, where lo^ is 
a high energy cut-off and a is a dimensionless coupling 
constant characterizing the strength of the dissipation. 
The two dimensionless couplings of the model are A/wc 
and a: in this paper we consider the region of parame- 
ter space, A/wc ^ a < 1, which includes the interesting 
case of a non-adiabatic bath. 

The Ohmic spin-boson model has been intensively in- 
vestigated over the last 10 years The main in- 
terest has been in understanding how the environment 
influences the dynamics of the two-level system and in 
particular how dissipation destroys quantum coherence 
[0,^1^-01 . Insight has been gained by exploiting the equiv- 
alence of the Ohmic two-state system to several other 
models, including the inverse square Ising model the 
anisotropic Kondo model (AKM) [^, and the resonant 



level model [p[-p^. The qualitative picture that has 
emerged can be summarized as follows. There is a renor- 
malized tunneling frequency A^, which depends on a and 
A, and which decreases monotonically with increasing a 
for fixed A. For a = 0, A^ = A (decoupled system plus 
bath) and for a > the renormalization of this scale 
increases dramatically as the dissipation strength is in- 
creased to 1: Ar/ci^c ^ (A/o^c)"'^^"'^^". Between a = 
and a — 1, there is a range of different behaviour from 
coherent oscillations at zero dissipation to damped oscil- 
lations at intermediate dissipation strengths and eventu- 
ally to incoherent relaxation at strong dissipation [|l|,D. 
At a > ac(A) k, 1 there is a "localization" transition at 
T — Q corresponding to a vanishing renormalized tunnel- 
ing frequency A^ = Here we address one aspect 
which has not been dealt with in a unified way in the 
literature, namely the meaning of universality and scal- 
ing in these models, in particular for physical properties, 
such as thermodynamic and dynamical quantities. We 
apply Wilson's numerical renormalization group (NRG) 
method to the AKM to calculate the specific heat, static 
susceptibility and dynamical susceptibility. The equiva- 
lence of the two models then allows us to discuss scaling 
and universality in the dissipative two-state system. 
The AKM @ is given by 
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where the first term represents non-interacting conduc- 
tion electrons and the second and third terms represent 
an exchange interaction between a localized spin 1/2 and 
the conduction electrons with strength J_l, J||. A local 
magnetic field, h, coupling only to the impurity spin in 
the Kondo model (the last term in Eq. ^ corresponds to 
a finite bias, e, in the spin-boson model. The correspon- 
dence between H and Hsb, established via bosonization 
[||, implies e = g^lBh, A. = p^j^ and a = (1 + ^f, 
where tan (5 = _Z£^^^ ^ jg phase shift for scatter- 
ing of electrons from a potential Jii/d and po = 1/21)0 is 



1 



the conduction electron density of states per spin at the 
Fermi level for a flat band of width 2Dq [Qj|J^. We choose 
LUc = 2-Do so that A — J± and measure all energies rela- 
tive to Dq = 1. Since we are interested in describing the 
Ohmic two-state system for A/lUc ^ a < 1, this requires 
in the AKM that poJ± < l-po-'^H for PqJ\\ < 1 (a ^ 1" 
case) and poJ± ^ l/{poJ\\)'^ for poJ\\ » 1 (a ^ case). 
The AKM for poJ± > poJ\\ will be dealt with elsewhere 
so in effect we consider only poJ± <C min{poJ\\, a). 

We solve the AKM using Wilson's NRG method 
In this procedure, a separation of energy scales is made 
by introducing a logarithmic mesh of k points fc„ = 
A~", A > 1, and transforming the Cka to a basis of Wan- 
nier states fna lH] at the impurity, with /oa = J2k ^ka, 

Ckfi, is tridiagonal in k-space. 
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such that He = J2k 

i-e. - E^E,?=oA-"/^(/I+i^/„M + /^-c.)- The 

Hamiltonian (0) in the new basis is now diagonalized 
by defining a sequence of finite size Hamiltonians Hn 
containing the first N Wannier states together with the 
impurity. One diagonalizes the rescaled Hamiltonians 

— N — 1 

H]si = A^~Hn which satisfy the recursion relation 
Hn+i = A^/^Hn + E^ifh+i^fN^ + h.c). This gives 
the excitations and eigenstates at a corresponding set 

N — 1 

of energy scales ujn defined by a; at = A ~ and al- 
lows the calculation of dynamic quantities at frequencies 
Lu ^ lon and thermodynamic quantities at temperatures 
ksTN ^ ^N- For example, the Fourier transform of 
x{t,T) = ^ie{t){[S,{t),SM]). is given by xi^.T) = 

X'(^,T) + ^x"(^,T) =^^E™,JM^.nPf^S^(^, 
where £,„ , are excitations of Hj^i , Ziq (T) the partition 
function of i/jv, and = {m\Sz\'n)N- 

Specifically, for the AKM we calculate: 
relaxation function S{lo) = _ i x 
rity specific heat C{T) = -d^hmpfd 
purity free energy is given by Fimp{T) = 



(a) the T = 

(b) the impu- 
where the im- 

'ksT In Z / Z() 



and Zq is the conduction electron partition function; 
(c) the local static susceptibility x'('^ = 0;^) corre- 
sponding to setting the g factor of the conduction elec- 
trons to zero. Under the equivalence, the operator 
of the spin-boson model translates to Sz in the 
AKM, so S{u}) gives the relaxation function for the spin- 
boson problem. We extract the local static susceptibility, 
Xsb{T) = -{l/p){d''\nZsb/de'),=o, of the Ohmic spin- 
boson model from x' {'^ = 0, T) at finite T. Finally, the 
decoupling of spin and charge degrees of freedom in the 
AKM allows identification of Fimp as the free energy 
of the spin-boson model (with bath contribution sub- 
tracted) and of C(T) as the corresponding specific heat. 

The energy level fiow diagram for some low lying 
rescaled energy states is shown in Fig. (]l|). We see that 
spin-rotational invariance, which is broken for Jx ^ J\\ 
at high energies, is restored below the low energy scale 
of the model, the Kondo temperature Tk{J^,J\\) @, 
leading to the well known isotropic strong coupling fixed 
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FIG. 1. The lowest rescaled energy levels for even A'^ it- 
erations for J|| = 0.443 and — 0.010 corresponding to 
A = 0.01 and a = 0.8 in the spin-boson model. The energy 
levels are labeled by the conserved quantum numbers, total 
pseudo-spin j ||l^ and total z-component of spin Sz- There 
is a crossover to the strong coupling fixed point at iteration 
Nc corresponding to A"''^''"^'''^ ~ = Tk- Spin rotational 
invariance is restored at low energies (e.g., the j = states 
with Sz = and Sz = ±1 become degenerate), so the states 
at the strong coupling fixed point can be labeled by total spin 
S as indicated. 



point at low energies (e.g. the lowest single particle 
states in Fig. (|l|), rji = 0.6555,7/2 — 1.976 agree with 
the A = 2 results of [^). A detailed analysis gives 
Tk ^ {A/ujc)^/^^~°'^ with A, a related to Jj_, J|| as above 
and Jj^ <C J|| and a prefactor also depending on a, i.e., 
Tk has the same dependence on a as the low energy 
scale Ar- The flow to the isotropic strong coupling flxed 
point holds for any initial anisotropy, corresponding to 
< a < 1, with the flow being universal for each a [[l5[ : 
the energy levels for flxed a and different A/ujc <C a 
may be shifted onto each other by a translation in N 
|]l6| , except for a small "transient" region TV '--^ 0-10 cor- 
responding to high energies w 3> near the cut-off uic- 
We see that the energy level flow is uniquely specifled 
by two parameters, A^ (equivalently Tk) which sets the 
crossover scale in Fig. (|l|), and the dissipation strength 
a (or equivalently the dimensionless initial coupling con- 
stant po^ii)- The universal flow of energy levels is the 
origin of the scaling, for flxed a and arbitrary A/ujc ^ a, 
in the thermodynamic quantities we discuss below jT^ . 

The inset to Fig. (|^) shows that the speciflc heat curves 
for a = 0.2 and several values of A all scale onto a univer- 
sal curve corresponding to a — 0.2. For different a one 
obtains distinct universal curves, C'a(T), Fig. (||). Scal- 
ing is valid for all temperatures in the range ksT <ti Dq, 
not only at low temperatures ksT ^ Ar — Tk- In Fig. 
(^, and throughout this paper, we scale the temperature 
by a/7 = SAr/iT^kj^- This follows from an exact result 
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FIG. 2. Universal specific heat curves, ° ^ ' , 

for the Ohmic two-state system for < a < 1. 
7 = liniT^o C(T)/T ~ ct/Ar is extracted from the fixed 
point analysis ||l5[. The symbols represent the temperatures 
at which the numerical second derivative of Fi^p has been 
calculated. The inset is for a — 0.2. 

for the Wilson ratio |H,|9| discussed below together with 
our definition Xsb(O) = l/2Aj.. The specific heat is linear 
in temperature for fc^T <^ A^. and < a < 1, with a lin- 
ear coefficient 7 in good agreement with values extracted 
from an analysis of the strong coupling fixed point . 
This is also expected from the Fermi liquid groundstate 
of the model. We note that the coefficient of the spe- 
cific heat changes sign close to q = 1/3, corresponding to 
the appearance, for weak dissipation, of damped oscilla- 
tions at frequency A^ ||,^. The peak in C{T)/T can be 
taken as a signature of a two-level system weakly cou- 
pled to bosonic excitations - it is absent for strong Ohmic 
dissipation. The strong dependence of C{T)/T on a for 
weak dissipation is also seen in other quantities, such as 
in the spin response 

The universal curves for the static susceptibility, 
Xa{T), parametrized by the dimensionless dissipation 
strength, a, are shown in Fig. (^. The inset shows that 
curves with the same a = 0.8 and different A scale onto 
the same curve for ksT <^ Dq. The thermodynamic cal- 
culation for Xa becomes inaccurate at low temperatures, 
ksT ^ Ar (discussed in detail in and we have 

to resort to an analysis about the strong coupling fixed 
point. This yields a finite susceptibility at T = which 
is accurate to within 1% for < a < 1 (see Table | below 
and fl^). At high temperatures, A^ ^ fc^r <C Dq, there 
is a dramatic difference in the approach of ksTxaiT) to 
its free spin value of 1/4 between the cases of weak and 
strong dissipation. In the former, the free spin value is 
reached very rapidly on increasing the temperature above 
Ar- The logarithmic terms characterizing the slow ap- 
proach of kgTxaiT) to the free spin value in the Kondo 



TABLE I. 7, Xakm and Rakm, extracted from the fixed 



point analysis. 


a 


A = Jx 


\ak77i 


7 


Rakm, 


10"^ 


10-" 


489.2 


3214.6 


2.003 


10"^ 


10-* 


4.78 


30.76 


2.04 


1/3 


0.01 


275.1 


1806.1 


2.005 


1/3 


0.1 


8.61 


56.4 


2.01 


0.7 


0.01 


1.90 X lO'' 


1.25 X 10^ 


2.002 


0.9 


0.1 


1.81 X 10" 


1.19 X 10^ 


2.003 



case are small for weak dissipation and only set in when 
a ^ 1^ . In this limit and for A/ujc ^ 1, the scaling 
functions for the Kondo problem are recovered for the 
specific heat and static susceptibility [|5|. 
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FIG. 3. Universal curves for the static susceptibility of 
the Ohmic two-state system, kBTxa{T), for < q < 1, with 
7 ~ a/Ar as in Fig. |. The inset is for a = 0.8 and A = 0.005 
(o), A = 0.05 (□) and A = 0.01 (o). 

A universal Wilson ratio for the AKM Q has been 
proven also for the Ohmic and non-Ohmic spin-boson 
models For the Ohmic case, the Wilson ratio for 

the spin-boson model, Rsi, — limT->o | (g^^^^^ '^'c'' ^ 2/a 
JlOt . This is related to the Wilson ratio, Rakm — 
limT^o|(^^^ - 2, for the AKM |] by Est = 

Rakm/ 01 Since "Xsb — Xakm 

/a lH and Xakm is the suscepti- 
bility of the AKM (with a g factor of 2 for the conduction 
electrons). Table | shows that Rakm — 2 (i?sb = 2 /a) is 
recovered for A /lOc *C a and < a < 1 . 

Dynamical quantities also show the universality dis- 
cussed above for thermodynamic quantities. The T = 
relaxation functions, S{uj), for dissipation strength < 
a < 1 have been given in ||]. A detailed analysis shows 
that these are universal functions oi lu/ Aj, parametrized 
by a: S = Sa{ijj/ A^). The case of a = 1/3, correspond- 
ing to the crossover between damped oscillations and in- 
coherent relaxation |||,^,^, is shown in Fig. (Q). Scaling 
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FIG. 4. Universal curve for the relaxation function Sa{i-o). 



in Sa{uj) extends to all frequencies uj <^ Dq and is not 
restricted to lu < A^. At high frequencies, uj 3> A^., we 
obtain S'„(w) - ^-(^-20)^ ^^ms x"{^) ^^(3-2^) 
Cs{t) = ([(T,(i),a,(0)]_) ~ 1 - ct2(i-") for D^^ « t « 
A"""^ with the a dependent exponents being accurate to 
within 0.1% for < a < 1. These results agree with 
short-time approximations Q and perturbative methods 
1^ in this limit. More importantly, they indicate, that 
scaling, in the sense discussed in this paper, can only 
be expected for fixed a. At low frequencies, lu <C A^, 
the Fermi liquid behaviour of x"{uj,T = 0) ^ aco, gives 
Csit) - -a/t^, with an a independent exponent for 
t» 1/A, 

Deviations from scaling, starting at high temperatures 
and frequencies, set in on increasing A/ujc for our finite 
bandwidth model. The scaling discussed here is valid for 
< a < 1 (0 < J|| < 00) as long as A/wc remains the 
smallest bare energy scale. 

The scaling and universality discussed above can be 
useful in interpreting experiments on dissipative two- 
state systems. The dissipation strength a can be deter- 
mined by fitting the data for some quantity, such as S{lu), 
to the appropriate universal scaling function ^'^(a;). The 
low energy scale, A^, can be extracted from the low en- 
ergy/temperature behaviour (e.g. from Sa{uj) by using 
the generahzed Shiba relation Sa{0) — 2a[xsbiT — 
0)]^ and Xsb{T = 0) = l/2A.r). In practice, this may be 
difficult for strong dissipation a > 1/3 since the scaling 
functions for different a will differ appreciably only for 
ksTjUJ > Ar- In this case an alternative is to extract a 



and Ar from 7 — t: k^a/3Ar and Xsb — l/2Ar. a and 
Ar can be more easily deduced for weakly dissipative 
systems a < 1/3, for which the scaling functions depend 
sensitively on a, even for /cgT, a; <C A^. 

To summarize, we have used the equivalence of the 
Ohmic spin-boson model to the AKM in order to study 



universality and scaling in these models. For anisotropics 
in the AKM, < J_l <C J|| < -l-oo, corresponding to 
dissipation strengths < a < 1 and bare tunneling fre- 
quencies A/lUc a in the Ohmic spin-boson model, 
the thermodynamic (dynamic) properties of these models 
are characterized by universal scaling functions of T/ A^ 
(to/Ar) which are distinct functions for different a and 
are independent of A — the latter entering only through 
Ar- As in the Kondo problem, the scaling functions are 
universal with deviations from scaling at high frequen- 
cies and temperatures arising from finite A/lUc- The dis- 
sipation strength in the Ohmic two-state system, just 
as the anisotropy in the AKM, determines the essential 
physics, in particular the renormalization of the low en- 
ergy scale Ar/ujc ^ {A/ujcy^^^^"^ and the form of the 
scaling functions. The perspective gained above may also 
be useful in understanding the highly anisotropic multi- 
channel Kondo models which arise in the context of single 
electron devices, and two-level systems in solids interact- 
ing with electrons . 
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